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We calculate the energy of the state closest to threshold for two and three identical, spinless 
particles confined to a cubic spatial volume with periodic boundary conditions and with zero total 
momentum in the finite-volume frame. The calculation is performed in relativistic quantum field 
theory with particles coupled via a Xcf)* interaction, and we work through order A®. The energy shifts 
begin at and we keep subleading terms proportional to IjL^, l/L® and 1/L®. These terms 

allow a non-trivial check of the results obtained from quantization conditions that hold for arbitrary 
interactions, namely that of Liischer for two particles and our recently developed formalism for 
three particles. We also compare to previously obtained results based on non-relativistic quantum 
mechanics. 


I. INTRODUCTION 

We have recently derived a quantization condition 
that, subject to some conditions, determines the energy 
spectrum for three particles in a cubic box [TJ [5] . The 
formalism is fully relativistic, and extends earlier work in 
non-relativistic effective field theories 011]. One applica¬ 
tion of this general result is to determine how the energy 
of the state lying closest to the three-particle threshold 
depends on the box size L. We have developed such a 
threshold expansion through 0{1/L^) in Ref. [S]. 

This expansion also allows a test of the formalism. 
Since the energy shifts are proportional to 1/T^, the 
kinematics becomes nonrelativistic for large enough L. 
Thus we can compare our results to those obtained us¬ 
ing non-relativistic quantum mechanics in Refs. EE]. 
This comparison is successful for the three leading terms 
(proportional to l/L" with n = 3,4,5), but turns out 
not to be useful for the 1/L® term. This is because the 
three-particle interaction enters at this order, and there 
is no a priori relation between a non-relativistic contact 
interaction and a relativistic three-particle amplitude at 
threshold^ Thus one can only use the 1/L® terms as a 
method for determining this relation, and not as a check 
on our threshold expansion. 

We have thus turned to perturbation theory (PT) as 
an alternative tool to provide a test of the 1/L® results 
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^ One indication that there cannot exist a simple a priori relation 
between the non-relativistic and relativistic three-body interac¬ 
tions is that both require regularization and are scale and scheme 
dependent, with different regularization schemes needed for the 
two theories. In fact, in the relativistic context the situation 
is even more complicated. This is because, as has been known 
for a long time (see, e.g. Refs. [HHU], and our recent discussion 
in Ref. ^), the standard three-particle scattering amplitude di¬ 
verges at threshold. To obtain a finite threshold amplitude one 
must perform subtractions, e.g. following the methods introduced 
in Refs. mis]. These, however, are not unique—many definitions 
are possible. 


from the threshold expansion. We pick the simplest inter¬ 
acting, perturbatively-renormalizable relativistic QFT— 
scalar \(j)^ theory—and determine the threshold energy 
shift through 0(A®), keeping terms scaling as 1/L" with 
n < 6 in the large volume expansion. Cubic order is 
sufficient to provide a non-trivial check on the thresh¬ 
old expansion developed in Ref. [5|. Furthermore, al¬ 
though the (j)'^ theory has no bare six-point vertex, there 
is an induced three-particle scattering amplitude start¬ 
ing at O(A^), and the subtraction methods developed 
in Refs. mm are tested at one-loop order by our cal¬ 
culation. In particular, as part of our calculation, we 
have worked out the subtracted three-particle amplitude 
at threshold through cubic order, using the subtraction 
defined in Ref. [5]. 

We have carried out the calculation both for two and 
three particles, so as to provide further cross-checks. In 
particular, we can compare the result of our perturbative 
threshold expansion for two particles with that obtained 
using the relativistic finite-volume two-particle formalism 
developed in Refs. [mils]. The agreement we find, as 
described below, gives us confidence in our methodology. 

We can also compare to the results for two parti¬ 
cles obtained using non-relativistic quantum mechanics 
in Ref. [Sj. Since the three-particle amplitude does not 
enter into this result, the comparison is unambiguous. 
We find a discrepancy in terms proportional to A^/L®, 
and make a suggestion for its source. 

The remainder of this article is organized as follows. In 
the following section we give an overview of our method¬ 
ology. Results for the two-particle threshold energy are 
then worked out in Sec. IIIII In Sec. Hyi we calculate the 
subtracted, three-particle scattering amplitude at thresh¬ 
old. In Sec. jVjwe determine the three-particle threshold 
energy and express it in terms of this subtracted ampli¬ 
tude as well as two-to-two scattering observables. We 
make the comparisons to prior results in Sec. |VI[ and 
present some conclusions. 

Three appendices contain technical details. Ap¬ 
pendix [^derives two identities for finite-volume momen¬ 
tum sums. Appendix |B| works out the s-wave scattering 
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length and effective range in \(j)^ theory. Finally, ap¬ 
pendix develops the threshold expansion for two par¬ 
ticles that follows from Liischer’s quantization condition 
out to 0{\/L^). 


II. OVERVIEW OF CALCULATION AND 
METHODOLOGY 

We consider the relativistic QFT described by the Eu¬ 
clidean Lagrangian density 

^ -h -f 

+ , (I) 

with (j) a scalar field. Although this theory contains no 
six-particle local interaction, a 3 —>■ 3 scattering ampli¬ 
tude M .3 is induced, as discussed below. This theory has 
a Z 2 symmetry, under which (/)—>■ —^, that forbids ampli¬ 
tudes involving an odd number of fields. This symmetry 
is also assumed in the general three-body formalism of 
Refs. [U [5]. We have included counterterms for wave- 
function and mass renormalization. These are tuned so 
that m is the physical, infinite-volume mass and that the 
residue of the infinite-volume propagator at the mass pole 
is unity. Note that we do not include an explicit countert¬ 
erm for the coupling Aq, preferring to work initially with 
the bare coupling and later describe its renormalization. 

To determine the energies of states close to threshold, 
we calculate even and odd particle-number correlators in 
finite volume. Choosing convenient overall factors, these 
are defined by 

C2{t) = (UrfMOr) , ( 2 ) 

C3(r) = {UrfMOf) , (3) 

where 

^p{t)= J d^x . (4) 

The subscript “L” indicates that the spatial volume is 
restricted to a cube of side L, with periodic boundary 
conditions applied to </>. Thus the allowed momenta are 
p = 21 : 71 /L, with n a vector of integers. We work in Eu¬ 
clidean space, with r the Euclidean time, which is taken 
to have infinite range. We choose to place all fields at 
zero spatial momentum since the threshold state in the 
absence of interactions consists of particles at rest. Thus 
our creation and annihilation operators will have large 
(0(1) in PT) overlap with the actual threshold state even 
in the presence of interactions. This is a convenience, but 
is not strictly necessary since all we need is for our oper¬ 
ators to have some overlap with the threshold state. 


The general form of these two correlators is known in 
terms of the eigenstates of the Hamiltonian of the theory. 
Assuming, as we do henceforth, that r > 0, we have 


C2{t) = 

\ ' 7 Ai52,TiT 

/ . ^ 2 ,r!,e ) 

nGeven 

(5) 

Csir) = 


(6) 


nGodd 


where 

= Ej^n - jm . (7) 

Due to the Z 2 symmetry, it is possible to separate states 
with even- and odd-particle quantum numbers into C 2 
and C 3 respectively. This implies that C 2 contains a 
contribution from the vacuum state, which [given the in¬ 
clusion of the exp(2mr) in its definition, Eq. (§] leads 
to a growing exponential of r. Similarly, C 3 contains an 
exponentially growing contribution from a single-particle 
state. Such growing exponentials might be problematic 
in a numerical simulation, but can be readily identified 
in an analytic calculation. 

Our aim is to pick out from the infinite sum of expo¬ 
nential contributions, those corresponding to the states 
nearest threshold. We know that there is only one such 
state for each correlation function, since there is only one 
such state in the free theory (with all particles at rest) 
and we are making an infinitesimal perturbation. For 
these threshold states, the quantities vanish as 

L —>■ 00 as a sum of powers of 1/L, up to possible loga¬ 
rithmic corrections—which in fact do not arise at the or¬ 
der we work—and exponentially suppressed terms. The 
latter, which behave as we neglect throughout. 

Such corrections are also dropped in Liischer’s general 
two-particle quantization condition and in our general 
three-particle formalism. As mentioned above, in this 
work we expand the threshold energy shift, AEj^„, in 
both Aq and 1/L, working through 0 (Aq) and 0 ( 1 /L^). 

As is well known (see, e.g.. Refs. [6l|7l Ha da), the 
leading contribution to the threshold energy shift is pro¬ 
portional to 1/L^. As we now explain, this implies that 
we can unambiguously pick out from the Cj (r) the contri¬ 
bution from the near-threshold state. Consider first the 
behavior of ^ from excited states. The lightest such 
state adds a minimal unit of relative momentum between 
two of the particles, so that AE 2 ,n = 2a;p — 2 m + 0{Xo), 
where ujp = + m^ with \p\ = 27r/L. For L ^ 1/m, 

the energy shift can be expanded as A£’ 2 ^„ oc oc 1/L^, 
and so we see that the excited states are parametrically 
separated from the near-threshold state, whose energy 
shift scales as 1/L^ as noted above. Thus there can be 
no avoided level crossings, and, in an analytic calcula¬ 
tion, we can unambiguously identify the excited state 
contributions. We stress that it is crucial to discard these 
exponentials before expanding in 1/L, so that the con¬ 
tributions do not become confused with the ground-state 
energy shift. Excited states involving more particles (e.g. 
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five-particle states in C3) are even more obviously sepa¬ 
rated, since then A.Ej^n ~ 2m, which does not vanish as 
L —>■ 00 . Similarly, far-subthreshold states, which con¬ 
tain less particles (e.g single particle states in C 3 ), have 
AEj^n ~ — 2 m, and the exponentials can also be easily 
separated. 

In light of the previous discussion, the method we use 
is as follows. We calculate C 2 and C 3 order by order in 
PT, and remove by hand the contributions from expo¬ 
nentially growing far-subthreshold states and from expo¬ 
nentially falling excited states. The resulting subtracted 
correlators we call Cj^thr- We know that these have the 
form 


Cj.thr = ^ythrC . (8) 

Thus if we expand in powers of r, and keep only the 
constant and linear terms, 

+ r[a,C,-thr(0)] +0(r2), (9) 

then the threshold energy shift is given by 




^rgj-,thr(0) 

C,.thr(0) 


( 10 ) 


The advantage of this method is that it allows us to keep 
track of powers of Aq in a straightforward manner. An al¬ 
ternative approach would be to identify the infinite set of 
perturbative diagrams leading to the exponential behav¬ 
ior in Eq. Q, with its associated renormalization factor 
This requires working to all orders in Aq in a sub¬ 
set of diagrams. We have used this alternate method as 
a check on our results, though we present no details here. 

The diagrams that we need to calculate to obtain 
Ai?,. thr up to third order in Aq are shown in Figs, [^and 
[^for C 2 and C 3 , respectively. The free propagators are 

{Mn)cl^{0,r2))x,=o = , ( 11 ) 


where, again, ujp = \/p‘^ + is the energy of a free par¬ 
ticle with a finite-volume momentum p. For each loop 
there will be a sum over spatial momentum restricted to 
the allowed finite-volume values. In addition there will be 
an integral over the Euclidean time of each intermediate 
vertex. For given values of each of the spatial momenta, 
this is simply an integral over exponentials, with the inte¬ 
grand depending on the time ordering due to the absolute 
values in the propagators, Eq. (11). Thus the integrals 
are trivial, but keeping track of all the time orderings 
is less so. We have used two independent Mathematica 
codes to ensure that all terms are included. After the 
vertex integrals, one can read off which terms are expo¬ 
nentially falling or growing, drop these by hand, and thus 
obtain Cj^thr(T). Expanding in powers of r leads to the 
results for C'j_thr(0) and 9TC'j,thr(0) [see Eq. ([^] and thus 
to AEj^thr [using Eq. ([To])]. 

The final stage of the calculation is, for each loop, to 
sum over the finite-volume momenta. This is done by 


> c 



(a) (b) (c) 

P 

>cx xocx 

(d) (e) (f) 



FIG. 1. Feynman diagrams contributing to the even particle- 
number correlation function. External particles have zero 
three-momentum. Diagram (i) also has a horizontally flipped 
partner, not shown. Examples of labels used in the text for 
loop momenta are shown. 


converting sums to integrals, which can be absorbed into 
infinite-volume loop contributions, plus a finite-volume 
residue leading to the desired power-law terms. The 
methods for converting sums to integrals are variants 
and extensions of those used in deriving the general 
finite-volume quantization conditions mnittsi and their 
threshold expansions [5]. The results we need are col¬ 
lected in Appendix [A] 

As a simple illustration of these methods, consider the 
diagrams of Figs. a) and[^a). By construction, both 
diagrams lead to C^“^(t) = 1. In particular, the contrac¬ 
tion factors cancel the 1/2 and 1/6 in Eqs. (H and (§, 
respectively. There are thus no exponentially growing or 
falling terms to remove by hand, and C'j_thr(T) = Cjir) = 
1 at this order. 


Since only the constant term in C/,thr is non-vanishing 
at leading order, the perturbative expansions of the quan¬ 
tities of interest can be written: 

00 

Q.thr(O) = 1 -(- ^ A^C];(L(0) , 

n—1 

(12) 

00 

a.Q.thr(0) = 5]A^[5.cj”L(0)], 

n—1 

(13) 

00 

n=l 

(14) 
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FIG. 2. Feynman diagrams contributing to the odd particle-number correlation function C 3 . External particles have vanishing 
spatial momentum. Not shown are diagrams obtained from (h), (j), (k) and (o) by applying the loop correction to the other 
vertex. An example of loop-momentum labeling is shown. Figs, (i)-(o) are also diagrams for the connected part of the three- 
particle scattering amplitude. We reference these diagrams, reinterpreted as infinite-volume scattering contributions, in our 
calculation of Sec. |IV| For that calculation the external particles are not assumed to be at rest. 


Inserting these expansions into Eq. (101, 


a41 = -5.C'(1(0), (15) 

= -drCflM + C(1(0) [d.C%iO)] , (16) 
AE% = + c%{0) [arCjy,,(o)] 

+ (cftl.(0) - t^i,l(0)') [5rC(l(0)] . (17) 


Thus a third-order calculation of AEj^thr requires a third 
order result for drCj^thiiO) but only a second-order result 
for C'j,thr(0). 

Another simple example of our methods is provided 
by the disconnected diagram of Fig. 0b). This diagram 
clearly has no dependence on r, so when we use Eq. (§ 
we find that (72 (r) oc Dropping this exponentially 

growing term leads to (72,thr = 0, so this diagram makes 
no contribution to the threshold energy. 

The final example we consider here is the connected 
diagram Fig. [^c), which involves a single vertex. All 
four propagators have vanishing spatial momenta, so we 
only need to do the integral over the vertex position, ti. 


The resulting contribution to the correlator is 
y{c) / ^ J- An 1 


cr'M = -i 


2 (2m)2 L3 t 


dri 1 


/ O pOD 

dTie*^^^+ / dri (18) 

-00 J T ^ 

(19) 


Sm'^L^ V 2m 


-b r 


Here there are no exponentially growing or falling terms 
to remove, and we find the first nontrivial result for the 
threshold energy shift 


^T^lthr(O) — ~ > 


( 20 ) 


as well as an 0(Ao) contribution to the constant term 


E^lr(O) = - 


1 

16m3L3 


( 21 ) 


We now describe three classes of diagram that we do 
not need to calculate explicitly, although they contribute 
to the Cj(T) at the order we work. This is because they 
either do not contribute to (7j,thr('r), or they lead only to 
changes in the overall normalization of the correlators. 
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FIG. 3. Examples of Feynman diagrams leading to mass and 
wavefunction renormalization in the correlator C 2 at the order 
we work. Similar diagrams involving the counterterms 5Z and 
SZm are not shown. 


Zj^thr, but not to the energy shifts AEj^thr- Examples of 
these diagrams are shown in Fig. [^for C 2 and Fig. |^for 
C3. 

The first class involves self-energy and counterterm in¬ 
sertions on the diagrams described above (i.e., those in 
Figs.[^and[^. All diagrams in Fig.|^and those of Fig. 
(a)-((^ are examples of this class. We first note that con¬ 
tributions involving tadpole diagrams, such as those in 
Figs. [§a)-(f) and Fig.ga), are cancelled identically by 
the SZm counterterm. This is because the loops are in¬ 
dependent of the external momenta and thus lead only 
to mass renormalization. The only subtlety is that SZm 
actually cancels the infinite-volume version of the tad¬ 
pole diagram, in which the loop is integrated rather than 
summed. However, since the difference between the sum 
and the integral is exponentially suppressed, scaling like 
g-mi, contribution remains in our 1/L expansion^ 

The other renormalization diagrams that contribute at 
the order we work are those that do not involve tadpoles. 
These are those shown in Figs.|^g)-(i) and Figs.|^b) and 
(c), together with related diagrams in which the renor¬ 
malizations appear on different external lines, and dia¬ 
grams containing the corresponding counterterms. Note 


^ The difference between loop sums and integrals is exponentially 
suppressed as long as there are no cuts through the loops in 
which, for threshold kinematics, all particles can go on shell. 
Several examples where this is not the case are described in sub¬ 
sequent sections. 


FIG. 4. Examples of Feynman diagrams whose contributions 
to C 3 lead only to renormalizations of the constant .^ 3 ,thr but 
not to changes in the threshold energy shift Ai? 3 _thr. Dia¬ 
grams similar to (a)-(c) but involving counterterms are not 
shown. 


that renormalizations appear only on external propaga¬ 
tors at this order. For all these diagrams the loop sums 
can be replaced by integrals—there are no power-law 
finite-volume residues. Thus, if the propagators were 
evaluated on mass shell, the contributions of these dia¬ 
grams would exactly cancel with those containing the ap¬ 
propriate mass and wavefunction counterterms (given our 
convention that the residue at the mass pole is unity). In 
fact, the external propagators, while evaluated at vanish¬ 
ing spatial momenta, are not Fourier transformed in time. 
Thus they contain excited state contributions. However, 
once these are removed, following our general prescrip¬ 
tion described above, we expect the general argument to 
hold and the contributions of these diagrams and those 
with the corresponding counterterms to cancel exactly. 
We have checked that this is the case by explicitly calcu¬ 
lating these diagrams. 

The second class of diagrams are those exemplified by 
Figs.gd), (e) and (f), as well as those obtained by hor¬ 
izontal reflection. The characteristic feature of this class 
is that there is a disconnected propagator joining two of 
the external fields at either the initial or final time, mul¬ 
tiplied by a “one-to-three” correlator. This factorization 
is maintained as higher order corrections are included. 
The disconnected propagator has no r dependence, and 
so provides only an overall factor. Thus, applying our 
methodology to this class of diagrams amounts to study¬ 
ing the three-particle threshold energy using the one-to- 
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three correlator. As noted above, this is a legitimate 
approach, since one can use any interpolating fields with 
the correct quantum numbers. In particular, this class 
of diagrams alone must give the same result for AAa^thr 
as that obtained from the full C 3 . Thus we can drop 
these one-to-three diagrams without changing the result 
for the energy shift. We have checked this argument by 
explicitly calculating all the diagrams in this class up to 
order Aq. Note that, since one vertex is needed to convert 
the initial single particle into three, only a second-order 
result for Ai? 3 _thr is obtained. 

The third and final class of diagrams are those exempli¬ 
fied by Figs.j^g), (h) and (i). Here one is effectively cal¬ 
culating the threshold energy shift using the “one-to-one” 
correlator, with the disconnected propagators at each end 
only changing the overall factor. Once again, this class of 
diagrams alone, analyzed using our method, must yield 
the correct result for Ai ?3 thrj and so can be dropped. 


III. TWO PARTICLE ENERGY SHIFT 


In this section we calculate the threshold energy shift 
for two particles. We work to 0(Ao) in PT and keep 
terms up to 0{1/L^) in the volume expansion. We have 
already obtained the contribution linear in Ao, Eq. (20), 


so we start here with the quadratic term, which arises 
from Figs. [i|d) and (e). 

For Fig. [^d), the form of the result depends on 
whether the loop momentum, p, vanishes or is non¬ 
vanishing. This is because one of the terms which enters 
the energy shift in the p = 0 case becomes an excited 
state exponential, to be discarded, in the case of non¬ 
vanishing loop momenta. The zero-momentum contribu¬ 
tion is 


C2 {t) D 


( V 

3 

3t 


8 m? 

4m 2 _ 


( 22 ) 


The term is the second term in the expansion of 


exp(—with AE^W^,. given in Eq. (20), and does 
not concern us here. For p ^ 0 we find the contribution 


C2(t) D 


^0 


1 ^ 1 


1^3 ^ 

p^O P 


- Sm^) 

2m{pPY 


+ T^ 
02 


-I- e 


— 2{ojp —m)T 


m 


(p2)2j ■ 


(23) 


The superscript A indicates that ultraviolet (UV) reg¬ 
ularization is required, although the choice of regulator 
is unimportant. The last term in the summand is from 
excited states, and is dropped in C' 2 ,thr(T), as explained 
in the previous section. The remaining two parts of the 
summand diverge for p —>■ 0 , and are converted to inte¬ 
grals plus finite-volume residues using the results in Ap¬ 
pendix]^ For example, the term proportional to r leads 


( 2 ^) 

to the following contribution to drC^ thr(O)' 




(0)D 


■E 


p#0 



64m^L3 


(24) 


(25) 


where we have used Eq. (A51 and introduced the short¬ 
hand J^= J d^p/{27rf. 

Turning to Fig. [^e), we note that p = 0 does not need 
to be treated separately. After dropping excited-state 
contributions, we obtain 


C'2,thr(T) D 


"^0 


A 

32to2l3 ^sE 
P 


-I- LOpin — 


( 2ma;4(, 


m) 




(26) 

Here the sum can be replaced by an integral, since the 
summand is nonsingular and we are dropping exponen¬ 
tially suppressed volume dependence. 


Combining the results from Eqs. (22), (23) and (26), 


and evaluating the sums using Eqs. (A5) and (A 8 ), we 


obtain 




drCiXM = 

A 2 X 

1 1 

5 

(27) 

8 m?-E? 256'iT'^m^L'^ 

256m^L® ’ 

cSLm = 

J A 2 



i0247r^m2L2 IQm^L^ 


(28) 


/(0.2) 

. -n A\ 


64:171^ 


+ 0 {L-^). 


Here A 9 is the one-loop integral defined in Appendix B, 
Eq. (B 6 ), while /(°’2) is the finite integral 


j(0,2) ^ 


f m 

1 

m 

Ipujp + m 

UJpjP 

a;4_ 


1 

Stt 


(29) 


Two comments are in order. First, the A 2 terms are 
exactly those needed to convert the factor of Aq multiply¬ 
ing the first-order results in Eqs. (20) and (21) into the 


renormalized A. The latter is defined in Eqs. (B7)-(B8), 
which we reproduce here for clarity: 


X = 32rmTa = Xo-A 2 XI +AsXl+OiX^), (30) 


where a is the two-particle scattering length. Second, we 
have truncated C>(I/L^), since higher order 

terms in 1 /L lead to contributions to of 0 ( 1 /L^). 

This is because C' 2 ^jhr( 0 ) multiplies 9rC';^\hr(0) 0 ^ l/L^ 

(S') 

when it contributes to as can be seen from 

Eq. ( [T^ . 

Now we move to third order. As is clear from 
Eq. ( |17| ), we only need to determine the term linear in t, 
order to obtain Furthermore, we 

can drop any contributions falling as I/T^ or faster. 

We begin with Fig. Bf). If both loop momenta vanish, 
the diagram is proportional to 1 /L® and can be dropped. 
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If only one loop momentum vanishes, we find (dropping 
higher order terms) 




1 


1 4to^ — 3^ 


lE 


1024to®L® ^—' (Mr 

p#0 

J 


ijpy 


40967r‘*m^L® 

1 


256to®L® 


/(1.3) 


(31) 


(32) 


where 


/(1.3) ^ 


jpLOp{Mp+m)p 




AMpfP' 


1 

271^ 



1 1 

f ^ V ^ 1 

5127772L® 1 


1 

(f ^ + 

I 1 

5127772^® 

\Jp WpP^ 

' 2 m®L® 


drC^AM D - 


819277 


This leaves Fig.[^i), and its horizontal reflection. Here 
we must treat the cases p = 0 and p 0 separately 
(see the label in the figure), since the separation between 
ground and excited states is different in the two cases. 
For p = 0, the contribution to the threshold correlator is 


drCi%M D 


1 


The UV divergent integral in Eq. (32) combines with that 
in the result from Fig. [^i) [given in Eq. p 8 | below] to 
give a term proportional to H 2 . This turns out to be 
exactly the contribution needed to convert the factor of 
Ag multiplying the r term in Eq. ( |2^ to A^. 

If both momenta are non-vanishing, the result factor¬ 
izes into a product of loop sum^ 


5.Cf,L(0) D I ^ V ^ 1 , (34) 

• (35) 


The only product of finite-volume residues that falls 
slowly enough to be included is 


(36) 


512to®L® 


1 


27rr 




4 — TT 

204877 ^ 771 ^^® 


m) 

rA 




512m®L® 07® 


(37) 

(38) 


The UV divergent integral leads to a renormalization of 
lower-order terms, as described above in the discussion 


(33) following Eq. (321. 


For p 7 ^ 0, the diagram contributes 

A 


drCS^M U - 


G(p, ^ = 


Ag 


1 


64m^L® L® 


E 

P¥= 0 ,q 


G{p,q) 






W 


OpU/qU^pg 


(1F2 - 


m?) ’ 


(39) 

(40) 


where W = (Mp+Lfg+iMpg and Wp, = y/ (p-l-^^ q- 77 ^ 2 ^ 
can replace the sum over q with an integral since the 
summand is regular, leading to 


A® 


^rC^Jhr(O) D g4^°2^3 ^3 E 

ppo 

■ G{p,q) , 1/(0) 




p2 


Aq 


jp 


Ait'll 


f(0) 

L® 


where 


/(p2)= Eg(p-',<7). 

J a 


(41) 

(42) 

(43) 


Terms involving a single hnite-volume residue multiplied 
by the integral give part of the contribution needed to 
convert the factor of A§ multiplying the and 1 /L® 

terms in Eq. (25) into A^. Terms involving two integrals 
contribute to two-loop renormalization, generating part 
of the A 3 A® term which converts Ag to A in the Ag/L® 
contribution to G 2 ,thr- 

The sums in both loops in Eig. [^g) and (h) can all be 
converted into integrals, and these integrals are exactly 
those obtained when the same diagrams are evaluated 
as contributions to infinite-volume scattering. It follows 
that these diagrams contribute only to renormalization 
of lower-order terms. 


® This factorization occurs because, in order to obtain a contri¬ 
bution proportional to r, the times of the vertices must satisfy 
0 < ri, T 2 , ra < T. The set-up is then essentially the same as 
in the calculation of a threshold scattering amplitude, for which 
we know, from using Feynman diagrams, that the contributions 
from the two loops factorize. This is true both in finite and 
infinite volume. 


Here we have assumed that the UV regulator maintains 
rotational invariance, and used Eq. (A51. The first term 
in Eq. (42) contributes to the two loop renormalization of 

^rG 2 \hr( 3 )> while the second completes the renormaliza¬ 
tion of the I term in Eq. (25). The third term completes 
the renormalization of the l/T® term in Eq. (251, leaving 
a finite residue which we now calculate. 

To determine f'{0) we must expand G for small p: 


G(p, q) = fo+p- qfi + P V 2 a + (P • ^26 + 


(44) 


where the fj are functions of q^. Performing the angular 
average over q, and picking out the term quadratic in p, 
we find 

f{0) = (^/2a(9^) + y/2b(9^)^ ■ (45) 

Carrying out the algebra, and the resulting finite inte¬ 
grals, we find 

i (ii - sie - /E) ■ <«> 







































The last term is part of the renormalization of the 1 /L® 
term in Eq. (25), as mentioned above. The first two terms 
complete the finite residues at the order we work. 


Adding the results from Eqs. (32), (36), (38) and (42) 
we find 


d fOl = ^ _ ^ 

T 2 ,thrl > 81927r4m4L5 40967 rTO 5 L 6 


(47) 


where the “renorm, parts” are the above-described con¬ 
tributions that convert Ag to A in lower-order terms. 

With all the results in hand, we can determine Ai? 2 ,thr 
through 0 (Aq/L®). We note that only the first and third 
terms in Eq. ([TtI ) contribute at this order, while both 
terms in Eq. ( |16p are needed. Writing the full result in 
terms of the renormalized coupling A we obtain 


_ A _ - J) 

Sm'^L^ 2bQTT'^m?L^ 81927r'^TO^L^ /^o'l 

QA2 >3 1 


IV. DIVERGENCE-FREE THREE-PARTICLE 
SCATTERING AMPLITUDE AT THRESHOLD 


The threshold energy shift in the three-particle case, 
depends, at C>(l/L®), on the three-particle scat¬ 
tering amplitude, M 3 . As is well know, however, M 3 is 
singular for certain choices of external momenta and, in 
particular, at threshold (see, e.g.. Refs. [SHII]). This is 
a well-understood physical singularity, described, for ex¬ 
ample in Ref. [T]. It implies that the I/L® term in AEg thr 
cannot depend on A43 itself, but rather on a subtracted 
version which is finite at threshold. In Ref. [1] we provide 
one possible definition for a subtracted amplitude, which 
we call the divergence-free amplitude. This definition is 
general, working for particles of any masses and both 
at and away from threshold. It is, however, a cumber¬ 
some definition to implement (e.g. involving an infinite 
number of subtractions, all but three of which are finite 
for degenerate particles at threshold). Our analysis of 
the threshold expansion of the three-particle quantiza¬ 
tion condition suggests instead using a simpler quantity, 
which we call A43^thr, whose definition involves only the 
minimum necessary subtractions This quantity is 
motivated and discussed at length in Ref. [5]. 

The purpose of this section is to calculate A43,thr to 
O(A^), so that we can express our perturbative result in 
terms of this infinite-volume quantity. Its definition is [5] 


A43,thr = lim 




M3^S — /o,5 — 

f I 


d^P „ 


Is ( 27 r) 

d^P 2 „ 


(P) 


(27r)3 Js (27r)3 


2 {P1,P2) 


(49) 


Here S is an IR cutoff, whose definition will be explained 
in the context of the following calculation. There are 
three subtraction terms, involving Jg, 5i and ^ 2 , respec¬ 
tively. Only the first two terms enter at the order we 
work, since ^2 = 0{X^). We give the definitions of the 
relevant parts of Iq^s and 5i below. 

The diagrams that contribute to M 3 at the order we 
work are those of Figs, [^i)-(o), now interpreted as Feyn¬ 
man diagrams in infinite volume. We first consider the 
0(Ag) diagram, Fig. [^i), which is also reproduced in 
Fig. ^a) along with the momentum labels we use. This 
gives a contribution that diverges at threshold, since the 
intermediate propagator goes on shell. The 5 prescription 
of Eq. (49) corresponds here simply to working away from 


threshold with general momenta, making the subtraction 
(here of /g), and then sending all external momenta to 
zero [ 5 ]. The contribution to the scattering amplitude is 


Mf = 


qZ — 777,2 _|_ 7g ’ 

q = (E - ujp - ujk,-p- k) 


(50) 

(51) 


where E = Up+uja+tOb, and we have used the vanishing 
of the total spatial momentum. The denominator of the 
propagator is 


-m^ = {E- W){E -W + 2uJup) 


(52) 


where ujkp = \J(p+kY + and, here, W = ujp+oJk-Mkp- 
This denominator vanishes when all the spatial momenta 
tend to zero, for then E and W both tend to 3m. It 
also vanishes for non-zero momenta li E = W (i.e. if 
uja + u}b = oJk + uJkp) and we must avoid such above¬ 
threshold divergent momentum configurations. The sub¬ 
traction term corresponding to this diagram is the lead¬ 
ing order part of Jg [S] 


Ai) 

° 2uJkpiE - W + ie) ' 


(53) 


By construction this has the same pole at if = IT as 
M^^\ so that the difference is finite (both at threshold 
and for above-threshold divergent momenta): 


M^ - /r = 


(i) _ 


A2 

^0 


2u)kp{E — IT -|- 2u}kp) 


(54) 


The final steps needed to obtain Algyhr are to sym¬ 
metrize over the external momentum assignments]^ and 
to take the threshold limit. Since the limiting value is 
independent of the choice of external momenta, sym- 
metrization gives a factor or 3 x 3, and the total con¬ 
tribution from this diagram is 


M 


9A2 


(7) __ 

3.thr 4^2 


(55) 


^ Since the particles are identical, interchanging a and b does not 
lead to a different assignment, so there are only three choices for 
each of the initial and final states 
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(a) 


(b) 



(c) 



FIG. 5. Feynman diagrams contributing to the three-particle 
scattering amplitude, M 3 . The momentum labels shown in 
(a) apply to all four diagrams. 


We next consider the diagrams in which the vertices in 
the tree diagram receive one-loop corrections. The full 
set of of these are those involving the s-, t- and u-channel 
loops shown in Figs. ib), (c) and (d), respectively, as 
well as the corresponding corrections to the right-hand 
vertex. The subtraction in this case is somewhat subtle 
so we provide a more detailed explanation. 

A key point in the following is that the intermediate 
propagator (with momentum q) is off shell, and so the 
loop corrections to the vertices at either end of this propa¬ 
gator differ in general from those which enter the on-shell 
scattering amplitude. This is true despite the fact that 
three of the four legs are on shell (e.g. those with mo¬ 
menta 6, c and k for the left-hand vertex). In fact, at 
one-loop order, the s-channel loops [e.g. Fig. [|:b)] do 
not depend on the off-shellness of g, while the t and u- 
channel loops [e.g. Figs, [^c) and (d)] do. 

This point is important because the general form of 
the subtraction term, Ig, replaces the loops at either end 
of the central propagator with on-shell scattering ampli¬ 
tudes, as well as changing the form of the propagator (as 
described above for the tree-level diagram). This feature 
of the subtraction is crucial, since it means that it is given 
in terms of physical quantities. The general form is quite 
complicated because it involves scattering amplitudes in 
all partial waves [5]. However, since we are interested in 
the subtracted amplitude at threshold, we need only the 
part of Iq that involves s-wave scattering, and we need 


this only close to threshold. Specifically, we have 


° 2 u}kp(E - W+le) 


(56) 


^)+0{b\X^), (57) 
6 m^ J 


where the superscript ^\ijk)” refers to the subfigures 
within Fig. [^ Here b* is the magnitude of the spatial 
momentum of b evaluated in the b + c CM frame, and 
b'* is the same quantity for 6'. The result (57) is sim¬ 


ply the threshold expansion of the two-particle scatter¬ 
ing amplitude, obtained using Eqs. (|B3|), (|B7|) and (|B14|) 


in Appendix [^ Note that this is expressed in terms of 
A = Ao(l — A 2 A 0 + ... ). Note also that the leading order 
part of Eq. ( [^ is the subtraction (53) that we used for 
the tree-level diagram. 

Evaluating Eigs. [^b), (c) and (d), together with the 
diagrams where the other vertex is corrected, and adding 
the result (50) from Eig.j^a) , we obtain 




M 2 (b* 


1927r2m2 ) 


X 


1 

+ ie 



1927r2m2 J 


(58) 


As noted above, only the t- and u-channel loops differ 
from the on-shell scattering amplitude. The part which 
differs is, close to threshold, proportional to t + u, where, 
for example, t = {b — k)^ and u = {c— fc)^ when the left- 
hand vertex is being corrected. Now t + u = +3m? — s, 

and s = (b + c)'^ is the same irrespective of whether g is 
on or off shell. Thus the difference between t + u on- 
and off-shell is equal to the difference between g^ on and 
off-shell, which is just g^ — m^. The key point is that 
this cancels the denominator of the propagator, leaving 
a finite residue in the threshold limit. 

We can now perform the subtraction, take the thresh¬ 
old limit, and multiply by the momentum permutation 
factor of 9, to obtain the final result from the diagrams 
of Fig. [5] 


M 


(ijk) _ 
3,thr 


9A^ 

Am? 


9A3 

967r2m^ 


+ 0 {X'^). 


(59) 


Note that we have written this result in terms of the 
renormalized coupling. 

The “bull’s head” diagram of Fig. il) also requires 
subtraction, in this case from the Si term in Eq. (49). 
At 0(Ag) the subtraction function is 


m{p) 


9Ag g(p)^ 
8 to {fPY ' 


(60) 


where H is the UV cutoff function discussed in Ap¬ 
pendix]^ Si is to be integrated over p with an IR cutoff 
IpI > 5, so that the integral is finite. The definition 
of the IR-regulated three-particle scattering amplitude. 
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A^ 3 ,i 5 is, in general, quite involved, as it requires apply¬ 
ing a cut-off both on loop momenta and also on the ex¬ 
ternal momenta. However, the external momenta scale 
as which implies that the IR cutoff they induce in 
the loop integral is weaker than that from the direct cut¬ 
off at Ip I = S. This means that we can set the external 
momenta to zero. Doing so, we find that 


M 


(0 _ 
3,5 “ 


9Ao f 3wp - 

16 Jp,5 ^3(^)2 


(61) 


where the subscript S on the integral indicates the IR 
cutoff. Combining with the integral of Si leads to 

an IR convergent integral in which one can set i5 = 0: 


iujp-m? 2H{p}'^ m? 
w3 m 

(62) 

We stress that this is a UV convergent integral that has 
a finite value, but that this value depends on the choice 
of cutoff function H. 

Next we consider the “twisted bull’s head” diagram of 
Fig. [^m). This diagram is convergent in the IR, so no 
subtraction is needed, and it can be evaluated directly at 
threshold. We find 


= - 


9Ag 

I6m2 


M‘’X = -. 6 A> 


(Yp 


3Ag 

IfiTT^m^ 


pi — rn^ + ie 


3 


(63) 

(64) 


The final diagrams involve intermediate single-particle 
states. These also have no IR divergences, require no 
subtractions, and can be evaluated directly at threshold. 
Since Fig.j^n) is a tree diagram, it is simple to evaluate, 
and yields 


M 


(n) _ 

3,thr 


A 


2 

0 




(65) 


The one-loop corrections to this diagram are given by 
Fig-d o) and the similar diagram where the right-hand 
vertex is corrected. Note that although this looks like an 
s-channel loop and might be expected to have t- and u- 
channel partners, in fact, no other diagrams exist. Eval¬ 
uating the loop we find that the sum of the two diagrams 
gives 


M 


(o) _ 

3,thr 


3Ag 1 

32 to 2 Jp ujp]? 


( 66 ) 


The UV divergence here is exactly that needed to convert 
Aq in (65) to A, up to a finite residue: 


■''^‘3,thr ~ 


A3 


8m? 327r2m2 


+ 0{XY. (67) 


where the last quantity is the contribution from 
Figs. [|;i)-(o), given by 


i^(Zmno) __^ A3 9A3C3 

3’**''' 8m2 8m2 ’ 


(69) 


with 


r m[iJ (^2 _ X] 

V (p2)2 


(70) 


We observe that, at the order we work, one could dispense 
with the cutoff function H, i.e. set H = 1 so that C 3 = 
0. This is appealing because it would remove the cutoff 
dependence from A 43 ^thr nt this order. However, this 
removal does not extend to higher orders in A. 


V. THREE PARTICLE ENERGY SHIFT 

In this section we determine the three-particle thresh¬ 
old energy shift, which (aside from renormalization) 
arises from the diagrams shown in Fig. The fully dis¬ 
connected diagram. Fig. [^a), was discussed in Sec. [lT| 
and leads to the leading “1” in Fq. (j^. The next 

seven diagrams. Figs. Sb) -(h), involve the interaction of 
only two particles, with the third being a spectator. Thus 
the results are the same as those for the corresponding 
two-particle diagrams, discussed in Sec. m except multi¬ 
plied by a factor of three for the number of pairs. Calling 
these contributions “disconnected”, we thus have, for the 
part linear in r, 

5rC;£L,disc(0) = 3 5.C'gL(0), n = 1,2,3,... (71) 

For the constant term, where we need only the terms lin¬ 
ear and quadratic in Aq, and only out to 1/^3 in the vol¬ 
ume expansion, it turns out that “connected” diagrams 
do not contribute. Thus the rescaled two-particle results 
are all that we need: 

(0) = 3 (b) > n = 1 and 2 . (72) 

Our remaining task is therefore to calculate the con¬ 
nected contributions to 9T-C'3_thr(0), which arise from 
Figs.[|;i)-(o). 

We begin with the tree diagram Fig. [^i). A simple 
calculation results ii{3 

SrcSLm ^ 3 ^ ■ (ri) 

One-loop corrections to this result are given by Figs. Hi) 
and (k) and their reflections. The contributions from van¬ 
ishing loop momenta scale as 1/L® and can be dropped. 


Combining results from all diagrams, we obtain 
9A2 3A3 


Al3,thr = 


4m2 327r2m2 


M 


(Imno) 

3,thr 


OiXY , (68) 


^ There is also contribution to proportional to l/T®, but 

this is beyond the order we need, as noted above. 
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Those from non-vanishing loop momenta give, after com¬ 
bining all four diagrams, 




L3 


256m^L® 
2 4 


SUk) 


(74) 


2m^ 


^;3 a;4(a;p-|-m) Wp(^)2_ 


(75) 


For the first three summands we can replace the sum 
with an integral at the order we work, while for the last 
we must use Eq. (A8). We also note that the first two 


summands are proportional to the integrand of A 2 [see 
Eq. (B6)], and in deed have the correct normalization to 


convert Aq in Eq. (It)^ to . Collecting all contributions 


and carrying out the resulting finite integrals, we find 

+ ^ + (76) 

Next we consider the bull’s head diagram, Eig. i:!)- 
Again the p = 0 contribution can be dropped. Eor non¬ 
zero loop momentum we find 




256m5L6 


m^(3Wp — TO^) 




(77) 

(78) 


Using Eq. (A8) we can write the sum as 





(3Wp - m^) 2 


m- 


+ 0 {L-^). 
(79) 

The integral appearing here is essentially the same as that 
in the subtracted bull’s head contribution to the scatter¬ 


ing amplitude, Eq. (62), except that the cutoff function 


H is absent. Combining results we find the total bull’s 
head contribution to be 


drCi%M 3 - 


377 


M 


(0 

3,thr 


3 C 3 


20487r4m4L5 


48Agm3L6 


128m3L® 

(80) 


The twisted bull’s head diagram of Eig. [^m) con¬ 
tributes 


o <-(3) —V — 

r 3.thrU J 128 to3L6 2,3/^ ^5 


1 


256Tr^m^L^ 


M 


(m) 

3,thr 


48AQm3L6 


(81) 

(82) 

(83) 


where in the second line we have converted the IR and 
UV finite sum to an integral and evaluated the integral, 
while in the third line we have used the result for the 


corresponding contribution to the scattering amplitude, 
Eq. (64). 

Finally we come to the diagrams involving single¬ 
particle intermediate states, Figs, [^n) and (o), together 
with the reflection of the latter. Since the presence of 
growing exponentials is a new feature of these diagrams, 
we give a little more detail here. The tree diagram of 
Fig. [^n) leads to 

^7 ^ ^0 f ^ ^2mr 13 rnT\ 

Dropping the enhanced exponential, the r term gives the 
contribution 




1 

384m^L^ 


M 


(n) 

3,thr 


48A§m3L6 ’ 


(85) 


where we have used the result Eq. ( |65| ) for the contribu¬ 
tion to the threshold amplitude from Fig. Sn). 

This by-now standard relation between the threshold 
amplitude and the 1/L® energy shift holds also for the 
one-loop diagram Fig. [^o). Including its reflection we 
find the contribution 




p^O 


( 86 ) 


Since the sum can be replaced by an integral up to 
0(1/L^) corrections, we see from Eq. ( |6^ that this con¬ 
tribution indeed equals A4 Shr/(48A3m3L6). 

The total result from all connected diagrams is 


Ag[^^C'Shr.co„„(0)]+A^[5.cg 


(3) / 

thr.conn ^ 


3A2 


3A3 


3A3 


32m^L^ 512'K'^m^L^ 20487rTO3L® 


M 


3A3C3 


128m3L® 48m^L^ 


(Imno) 

3,thr 


+ 0(A"), (87) 


where is given in Eq. (69). We observe that 


the XqJ/L^ terms cancel between Eigs.|2y), (k) and (1). 
Combining this with the results of Eqs. (|71|) and (72) for 


the disconnected diagrams, and using the relation (68) 
between and together with expressions 

(15 )-(|l7|) for the energy shift, we obtain 


3A 


3X^1 


8 m?L'^ 256'iT'^m^L'^ 


3A3 


3X^(1^+ J) 9A2 
81927r"‘m"‘L3 256m^L^ 256Tr'^m^L^ 

3X^0^ A43_thr 


128m3L® 48m3L® 


+ 0 {XyL^,l/L^)- (88) 


This is the main result of this article. We recall that —A 
is defined to equal the threshold two-to-two scattering 
amplitude, and is related to the scattering length via A = 
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32m7ro. Finally we observe that the dependence on the 
cutoff function H cancels between the last two terms, as 
must be the case since the finite-volume energy shift is a 
physical quantity. 


VI. COMPARISONS AND CONCLUSIONS 


To compare our results for the threshold energy shifts 
to those in the literature, it is convenient to re-express 
them in terms of the scattering length. For the two- 
particle case we obtain 


af;. 


2.thr — 


47ra 

rriL^ 


1 - 


al 


ttL 


7r2F2 


(6) 

-^+0{L-^) 




ina 


Stto 

9 


32a2 

3m 


C>(a^ 


(89) 

(90) 


The 1/F^, 1/i^ and 1/L® contributions agree with those 
obtained in Refs. [HI [H] (Reference [7] did not consider 
this quantity.) To aid comparison between the 1/L® 
terms, we rewrite in terms of the effective range r 


given in Eq. (B151: 


,( 6 ) 


a, = —- 


-h 


+ 0 {a^). 


(91) 


The latter form agrees with that obtained in Appendix [C| 
from expanding Liischer’s quantization condition m to 
0{1/L^). It is, however, in disagreement with 


,( 6 ) 


(Ref. [5]) = 


S-TT^o^r 


+ 0 {a^) . 


(92) 


The disagreement, which is proportional to X'^jL^ in 
Ai? 2 ,thr, appears to arise from the use of a non-relativistic 
dispersion relation at the last stage of the calculation, as 
discussed in Appendix We stress that a^r starts at 
O(X^) in PT, so that the result from Ref. [B] does contain 
a term scaling as /LX' , but with a different coefficient 
from that which we find here. 

We note that Liischer’s quantization condition holds 
for general scalar held theories, including effective held 
theories with arbitrary higher-order couplings. Thus the 
form of that holds in all such theories is that given 
in Eq. (91), i.e. i n te rms of a and r. By contrast, the 
form given in Eq. (90) holds only for A(/)^ theory. 


Our result for the three-particle energy shift is 


AE: 


'3,thr — 


127ra 

rriL^ 


1 -- 


aX 

ttL 


7r2L2 ) 


where 
(6) 127ra 

tto = - 

m 

127ra 


37ra 
T ’ 


32a^\ reSa^Tr^Ca M 


m 


m 

ina 

m? 




+ - 


m^ 

reSa^TT^Ca 

m? 


(6) 

(93) 

(94) 

(95) 


3,thr 


48mX 

AIa,thr 
' 48m3 


As for the two particle case, the second expression holds 
for general interactions whereas that in terms of a alone 
is special to Xtj)'^ theory. The 1/F^, 1/L'* and 1/L^ con¬ 
tributions agree with those obtained in Refs. HE]. As 


for 


,( 6 ) 


, Ref. H finds 


247r^a^r 


if ^(Ref. H) = -h r] 3 {fi) + 0{a'^), (96) 


while Ref. El quotes 


.(6) 


(Ref. [7]) = 




+ D + Oia^). (97) 


Here rj^ (p) is a three-particle contact interaction while D 
is the “three-body scattering hypervolume”. Both char¬ 
acterize a local three-particle interaction within the non- 
relativistic context of their respective calculations. The 
scale ^ in 773 is a renormalization scale, and any depen¬ 
dence on this scale must cancel out in the energy shift. 
The results of Ref. [Bj show, however, that this depen¬ 
dence enters only at 0 {a'^). 

While 773 and D are non-relativistic analogs of our 
threshold amplitude A 43 _thr, there could, in general, be 
finite differences between these quantities. Indeed equat¬ 


ing the three results in Eqs. (95), (96) and (97) gives 


relations between these quantities. However, given that 
our relativistic calculation gives a different result for 
a' 2 '^ from that obtained using non-relativistic methods in 
Ref. [Bjj it is not clear whether the non-relativistic three 
particle results apply at 0[1/L^) in a relativistic theory. 

Irrespective of these considerations, our result allows 
us to check the threshold expansion that we obtain in 
Ref. [5| from our relativistic three-particle quantization 
condition HE]. We find complete agreement, giving us 
further conhdence in the quantization condition. Fur¬ 
thermore, the perturbative calculation carried out here 
has provided a first explicit verihcation of the details of 
the subtractions needed to define a finite three-particle 
scattering amplitude at threshold. 

It would be interesting to push this calculation to one 
higher order in A, so as to allow a check of the A'^/L® 
terms that arise in the threshold expansion of Ref. [5] . 
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Appendix A: Identities for finite-volume sums 

We collect here various results needed in the main 
text to convert sums into integrals plus a finite-volume 
residue. Throughout we drop terms that are suppressed 
exponentially, i.e. as exp(— tttL)/L” for n > 0. We make 
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use of a cutoff function H{p) that was introduced in 
Ref. [1] 1^ The relevant properties of H are that it equals 
unity for = 0, vanishes for > 16 to^/ 9, and inter¬ 
polates smoothly in between. In addition, all derivatives 
of at p^ = 0 and at = 16m^/9 vanish. 

We first consider the generic sum 


Si 







(Al) 


where I is a geometrical factor, equal to Zoo(l,0) of 
Ref. [12]. Combining these results we obtain 


Si 


fif) T/(0) 


fjo) 

L3 




(A5) 


In other words, the sum can be replaced by an integral 
with exponential accuracy aside from a residue consisting 
of a 1/L and a l/L^ term. 

The second sum we need is the generic “l/(p^)^” sum 


As always in this appendix, the sum is over finite-volume 
momenta, in this case excluding p = 0 where the sum¬ 
mand diverges. The function f{p'^) is assumed regular 
at p = 0, so that the corresponding integral is IR conver¬ 
gent in three dimensions. In general, however, the sum 
is UV divergent, and must be regularized in some way, 
as indicated by the superscript A on the sum. All we 
need to know about this regularization is that it involves 
a cutoff scale A 4 to/3, so that H = 0 at the cutoff 
scale. 

To proceed, we rewrite the sum a^ 


^7 



mHW\ 

) 


/(O)g(p)^ 


Jp 


(A2) 


Q 1 V g(P^) 

^ L 3 ^ (^)2 ■ 

p^O ’ 


This sum can be rewritten as 

L3 


a _ 1 Y^ g(p)~g(o) , g(o) 

“ 7,3 (p2 V2 7,3 

p^O 


L 3 ^ (^)2 ■ 

p#0 '' 


(A6) 


(A7) 


The first term now diverges only as 1/p^ in the IR, so 
that, using the result ( |A5[ ), the sum can be replaced 
by an integral up to 0{i/L) corrections. The sec¬ 
ond term is proportional to a sum over integer vectors 
J = which equals Zoo{ 2 , 0 ) in the notation 

of Ref. [121 . Thus we obtain 




g(^)-g(o) 
Jp {P^r 


Lg{0)J 

IGtt^ 


+ 0{1/L). 


(AS) 


Appendix B: Scattering length and effective range 


Given the properties of H, the combined summand on the 
first line is regular at p = 0. Thus, if one adds the p = 0 
term to the sum, it can be replaced by an integral up to 
exponentially suppressed corrections. Thus the first line 
becomes 


_ fiO)H{p}^ 

Jp \ P^ P^ 


m 

L3 


+ 0{e-"^^). (A3) 


In the sum-integral difference on the second line of 
Eq. (A2) the cut-off A can be dropped since the 
summand-integrand is regulated in the UV by H^. This 
difference is then nothing other than a regulated I/p^ 
sum. Using a result derived in Ref. [5|, the second line 
can be written 


J/(0) 

dTT^L 


+ 0(e-™^). 


(A4) 


In the main text we need the perturbative expressions 
for the s-wave scattering length a and the effective range r 
in the scalar X(j)‘^ theory. The standard definition of these 
quantities in terms of the s-wave phase shift <5o (using the 
nuclear physics sign convention for a) is 


tan do (g) 


= —a I 


rag 


+ 0(q^) 


(Bl) 


Here q is the magnitude of the spatial momentum of each 
particle in the CM (center of mass) frame, which is re¬ 
lated to the total CM energy by E 2 = 2-^q^ + m^. The 
phase shift itself is related to the two-particle s-wave scat¬ 
tering K-matrix (defined as the angular average of the full 
K-matrix) by 

A2.s(g) = 167rU2 . (B2) 

q 


® H{p) is actually a function of but we use p as the argument 
for brevity. We note that the derivation can be made without 
introducing H, but doing so allows us to borrow a result from 
Ref. [J. 

^ The appearance of in the extra terms is simply to match 
the results that appear in the threshold expansion of the three- 
particle quantization condition For the derivation we could 
equally well use H alone. 


The K-matrix is related to the usual scattering amplitude 

by 


1 ^ 1 _ . q 

M2,s{q) A2,s(g) *167r£’2 ■ 

It follows that the scattering length is given by 


(B3) 


(B4) 


327rma = —/C2,s(0) = —A42,s(0). 
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Through one-loop order, evaluating the diagrams of 
Fig-Qc), (d), (e) and the u-channel version of (e) [treat¬ 
ing these as Feynman diagrams for infinite-volume scat¬ 
tering] one hnds 


Expanding in powers of we obtair 


PV 


1 




Jp 


(B12) 


Al2,.(0) = -Ao{1-AoA2+0(A3)} , 


fA 

■ 1 1 ■ 

h 

AujpjP 2^3 _ 


(B5) 

(B 6 ) 


Here we have done the integral over the loop energy vari¬ 
able po so as to write the result in the same form as those 
we obtain in the main text. The first term in square 
braces arises from the s-channel loop, while the second 
is the sum of those from the t- and u-channel loops. As 
usual, the superscript A indicates an (unspecihed) UV 
regulator. 

We hnd it convenient to adopt a physical renormaliza¬ 
tion condition for the coupling constant, defining it in 
terms of the scattering length as 


A = —A42,s( 0) = 327rma (B7) 

= Aq (l — A 0 A 2 -I- AqA 3 -|- C>(Aq)) , (B 8 ) 


where A 3 is the third-order coefficient that we do not 
need explicitly. 


To obtain the effective range we expand the K-matrix 
away from threshold in powers of q^. Specifically, using 
Eqs. (Bl), (B2) and (B 8 ) we find 


The first term reproduces the s-channel contribution to 
A 2 [cf. Eqs. (B5) and (B 6 )]. The second term leads to 
the following contribution to the derivative: 


,dIC 


2,s 


dq^ 


D - 


13^=0 


IGtt^ 


Combining Eqs. (|B10|) and (|B13|) we find 
2 


dq^ 


9^=0 


127r2 


o{K) ■ 


(B13) 


(B14) 


Comparing to the definition of r in Eq. (B9) we obtain 
our final result 


rm = + X— + C>(A). (B15) 

A 37r 

We stress again that the strange looking 1/A term is of 
purely kinematic origin. The non-trivial result of the 
one-loop calculation is the constant term. 


Appendix C: 1/L® term from Liischer’s quantization 
condition 


m 


,dJC^ 

dq^ 


= -A 

1 , rm 

^ F A- 

q ^—0 

2 647r 


(B9) 


where the first term in the square braces is of kinematic 
origin, coming from the expansion of E 2 in Eq. (B2 1 . 

Dependence on q^ first arises at one-loop order. For the 
sum of t- and u-channel loop-diagrams, a straightforward 
calculation finds that the scattering is pure s-wave, and 
yields the contribution 


Liischer’s original work on the two-particle thresh¬ 
old energy shift presented explicit results only up to 
0{1/L^) [12 ■ To compare to our perturbative result we 
need also the 1 /L® term from the general quantization 
condition of Ref. m- To our knowledge, the explicit re¬ 
sult for this term has not been presented elsewhere, so 
we determine it here. 

We start from the quantization condition in the form 
given in Ref. [T], itself adapted from the form derived in 
Ref. |14j . For vanishing total momentum, the condition 
is [Eq. (96) of Ref. [T] with fc = 0]: 


,dJC2 


m 


dq^ 


D - 




487r2 


(BIO) 


Note that the IC 2 = M 2 for these diagrams, since there 
are no physical cuts. 

The s-channel loop diagram does have a cut, so IC 2 
and A 42 differ. We elect to calculate the former, which 
requires using the principle value prescription for the pole 
that remains after doing the energy-component integral. 
The contribution is purely s-wave, and gives 


A2 f 

/C 2 (s-channel) = -^PV / 

8 Jv 


ujpijr^ - g2) ■ 


(Bll) 


det (/C/-1 -p F) = 0. (Cl) 

Here IC 2 and the kinematical function F are matrices 
in angular-momentum space. IC 2 is diagonal, whereas 
F contains off-diagonal elements. Close to threshold, s- 
wave scattering dominates, and one can show that it is 
sufficient to truncate the quantization condition to the s- 
wave alone, up to 0 ( 1 /in the energy shift, at which 


One must evaluate the principle value integral for > 0 to 
obtain this result. To obtain the correct result with < 0, 
one must use a modified PV prescription that yields analytic 
dependence on q^ [I]. 
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point the £ = 4 amplitude is needed. Thus, for our pur¬ 
poses, the condition reduces to 


truncating here at the a® term. It follows from Eq. ( |C6[ ) 
that has a similar expansion 


JC:;l+F, = o, 


with 


F, = 






_ 8uj^{E 2 - 2ujp) ' 


(C2) 


(C3) 


Here H is the cutoff function discussed in Appendix [Al 
E 2 is the two-particle CM energy, and the PV pole pre¬ 
scription is a generalized principle-value prescription in¬ 
troduced in Ref. [T] . We stress that this condition is iden¬ 
tical to the truncated Liischer quantization condition of 
Ref. [13] up to exponentially suppressed terms. 

In our companion analysis of the three-particle quanti¬ 
zation condition near threshold, we derive the threshold 
expansion of Eg (called 2mFoo in that work) |5|: 


Eg = 


AE 2 


1 


q'^L^ Att^L 
{q^L^fK. 

(47r2L)3 


q^L^J 

■ (47r2L)2 

0(L-3)}, 


(C4) 


where q is dehned as in Appendix the geometric 
quantities I and ££ are described in Appendix and 
JC = Foo(3,0) is a third such quantity (evaluated in 
Refs. [31 HD). The other result we need to apply the 
quantization condition is the expansion of 1 /A 2 ,s, which, 
from Eqs. (iBlj) and (B2|, is 


’<^2,1=-- 


1 




167raF2 . 

The kinematic relation we need is 

q^ = mAE2 + AF 2 /4 , 


(C5) 


(C6) 


2 

n—3 


(C8) 


Inserting the expansions of Eg and IC 2 ^ into the quan- 


tization condition (C2), and using Eq. (C81, we find 


5 ( 3 ) = 47]-a ^ 
^ _ _a 

6(3) 7r 

W) 

6 (6) 


2Trra^+(^^^ (-1^ + 3IJ -/C). 


(C9) 

(CIO) 

(CIl) 

(CI2) 


Converting this to a result for AE using Eq. (C6) we 
obtain the desired results 


t ^ 6(”) 

(n) _ e 


= - (n = 3 - 5), 

m 


(6) 

,(3) 


W) 


ira 

vn? 


(CI3) 

(CI4) 


“■2 arises from the expansion of the rel- 


The last term in a 
ativistic form of the energy, i.e. from the second term in 
Eq. (C6). The results for o„ agree with those of Ref. [5] 
aside from this term, suggesting that the source of the 
discrepancy is the use of a nonrelativistic dispersion re¬ 
lation in that work. 

In the main text we work only to cubic order in A. At 
this order we have 


where AEo = En — 2m. We are interested in the solution 


to Eq. (C2) near threshold, so we expand AE as 


AE 2 = y ^ 

Tn 


n—3 


(n) 


4 ‘’ = — 

m V 


OiX'^ 


(C7) where a and ra^ both start at 0(A). 


(CI5) 
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